Intrinsic non-linear conduction in the super-insulating state of thin TiN films 
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We investigate experimentally the electric transport at the insulating side of the superconductor 
to insulator transition in thin TiN-films. At temperatures T > 50 mK we observe an Arrhenius- 
type conductance, with an activation energy depending logarithmically on the sample size. At high 
bias the current voltage (I-V) characteristics display a large current jump into an electron heating 
dominated regime. For the largest samples, and below 50 mK we observe a low-bias power law 
/ oc V a characteristics with an exponent a > 1 rapidly growing with decreasing temperature, which 
is expected for a binding-unbinding crossover of the charge-Berezinskii-Kosterlitz-Thouless type. 

PACS numbers: 74.25. fc, 73.50.-h, 74.25.Dw. 74.78.-w 
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An intriguing low-temperature highly resistive state 
forms in disordered InO and TiN films approaching the 
disorder-driven superconductor-insulator transition (D- 
SIT) from its insulating side [HQ. A characteristic fea- 
ture of this state is the switching character of the current- 
voltage (I-V) characteristics: at certain magnetic field 
( independent threshold voltage Vhl the current jumps 
up over several orders of magnitude. It was further found 
that this highly resistive state in TiN can show hyperac- 



tivated temperature behavior of the resistance 

Comparing the data on the TiN films and the two- 
dimensional Josephson junction arrays (JJA) [1 . and 
exploiting the duality between the supcrconductingand 
insulating states in JJA 0, it was proposed 0, IH 
that in disordered TiN films a distinct superinsulating 
state forms, which is a low temperature phase of the 
charge Berezinskii-Kosterlitz-Thouless (BKT) transition. 
Within the JJA-model conduction in the insulating state 
is caused by tunneling of Cooper pairs from one island to 
the next, leaving an opposite charge behind. At low tem- 
peratures long range Coulomb interactions start to dom- 
inate the conduction process: a single chargc/anti-charge 
pair generates polarization charges all over the sample, 
which requires an energy that increases in two dimensions 
logarithmically with the size of the sample (provided the 
electrostatic screening length exceeds the sample size). 
Hence the thermal activation energy becomes dependent 
on the sample size 0, H| ■ The tunneling requires a cas- 
cade energy relaxation via the emission of electric dipolc 
excitations and phonons. The microscopic mechanism 
of superinsulation was proposed to be the suppression of 
this relaxation and thus blockade of the tunneling current 
below the charge BKT transition Q . 

A detailed study of the current jumps in InO was un- 
dertaken in [l(|. It was found that the I-V isotherms 
measured on the magnetic field-tuned insulator at B = 
1 1 T above the B value of the maximum of the magne- 
toresistance peak that, for this sample, is at B = 8T, 



and slightly exceeding the upper critical field, B C 2, are 
described excellently by the standard heating instability 
model 11 -ill- The latter was adapted [l4[ to strongly 
disordered materials by substituting into general expres- 
sions of [ll|-[l3j the electron-phonon relaxation rates 
calculated for strongly disordered conductors in [iUH 
(see [l2| for a comprehensive review). 

The findings of [Io| presented convincing evidence that 
at high fields Cooper pairing the InO films is destroyed 
and these films behave like conventional strong insula- 
tors. On the other hand, fundamental questions calling 
for detailed experimental investigation remain: (i) The 
nature of the high-resistive state near the D-SIT where a 
Cooper pair insulator (CPI) forms, (ii) The role, if any, 
of the long range Coulomb interaction in formation CPI 
properties, and (iii) The role of charge BKT in formation 
of the high resistive state of CPI. 

In this Letter, we address these questions. We measure 
the activation energy in the CPI-statc, and find it indeed 
depending logarithmically upon the sample size. Our 
data reveal a charge BKT-like cross-over, which marks 
the formation of a low temperature state with highly non- 
linear I-V curves. Below 50 mK I(V) follows a power-law 
/ oc V a behavior with an exponent a rapidly growing 
with the decrease of temperature, which provides evi- 
dence for super-insulating behavior at very low tempera- 
tures. Finally, we uncover a broad distribution of switch- 
ing voltages from the HR- to the low-resistive (LR) state 
at low temperatures in a strong contrast to switching 
governed by the overheating instability observed in [l(| • 

We measured the /-^-characteristics on the 5 nm thin 
disordered TiN films on the verge of the D-SIT similar by 
their parameters, to those studied in 0, 01, El ■ The data 
presented below are taken on two samples: the sample S 
which is superconducting at zero magnetic field but can 
be driven into an insulating state by applying a small 
magnetic field, and the sample I, remaining insulating 
even at B = 0. Magnetic field was taken perpendicu- 
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FIG. 1. (color online) (a) Ra vs. 1/T for Sample S at different 
magnetic fields B. Inset: Activation energy To vs. B of sample 
S. Line is guide for an eye. (b) Sample I for different sample 
lengths at zero magnetic field. Insets: To vs. sample size L 
of Sample I at zero magnetic field (top) and optical image of 
the 5x5 /zm device of sample I (bottom). Solid lines on both 
main panels are fits to Arrhenius law. 
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FIG. 2. (color online) I-V characteristics (a) of sample S 
at B = 0.9 T and T = 100,80,70,60,50,45,42,40,37 and 
32 mK from top to bottom; and (b) of sample I at B = T 
and T = 80,70,60,50,40,32 and 24 mK and L = 240 /mi. 
The voltage is swept from negative to positive values, and 
only a positive voltages are presented. Solid lines indicate 
the slopes corresponding to different values of power a in the 
I oc V a dependence. Insets in (a) and (b) display a(T) for 
both samples and I-V in linear scale for Sample I, respectively. 



lar to the sample plane. The sample S has a Hall-bar 
form. All I-Vs were measured on a stripe connecting the 
two most distant contacts (13 squares corresponding to 
650 /im and 50 /im in length and width, respectively). 
The transition temperature at zero magnetic field is be- 
low 100 mK. On a different chip of the same wafer, the 
sample I was patterned into the six independent squares 
of different side lengths, L (L = 0.5, 2, 5, 20, 240, and 
500 yum), and attached to two wide and thick gold pads, 
to provide a uniform current distribution and low contact 
resistances (see inset in Fig. lb). After the patterning 
process these samples were found insulating for T — > K 
at B = 0T. The results presented below are based on 
the detailed analysis of about 4000 dc I-V characteris- 
tics measured by the two-probe technique. Bias voltages 
were provided by a dc voltage source Yokogawa 7651 via 
a voltage divider (1000:1), with a highest sweeping rate 
of 1 mV/min. Currents were measured with DL1211 and 
Femto DDPCA-S current preamplifiers. In the supercon- 
ducting regime current biased four-probe measurements 
were carried out; in the insulating regime two-probe volt- 
age biased measurements were sufficient as the contact 
resistance was negligible below 1 K. The used setup en- 
abled us to measure resistances exceeding 100 Gfi. 

We first discuss the effects of a magnetic field on the 
sheet resistance of the S-sample. In the Fig. [lji we 
show log Rq vs. 1/T plots where each data point de- 
picts the zero-bias resistance extracted from the dc I- 
y-characteristics in the linear response regime. The 
superconducting state of the Sample S is very fragile 
and suppressed by small magnetic fields B < 5 mT. 
Solid lines represent fits to an Arrhenius law (i?rj(T) = 



Ro exp (Tq/T), where IcbTq is an activation energy). At 
the lowest temperatures R(T) deviates from the Arrhe- 
nius behavior and for T < 50 mK no linear resistance 
could be extracted from the i-V^-curves (see below). In 
parallel magnetic field qualitatively similar behavior was 
found with two orders of magnitude smaller typical resis- 
tances (not shown). The magnetic field dependence of T 
exhibits a nonmonotonic T-dcpendcnce with a maximum 
around 0.9 T similar to earlier experiments [1, [H HH , and 
is plotted in the inset of Fig. QJi. Employing Joscphson- 
junctions network (JJN) model Q to Tq(B), one can es- 
timate ds = 28 nm for the size of the Joscphson junctions 
plaqucttes for Sample S . 

In Fig. [TJd we present logi?rj(l/T) of squares with in- 
creasing side length L for Sample I. The i?rj(l/T) curves 
display again Arrhenius behavior at higher temperatures, 
but at low T striking differences occur: while the largest 
sample appear to follow the Arrhenius law down to the 
lowest temperatures, for which a linear resistance can be 
extracted, the smaller samples exhibit a saturation of Rjj 
at low T, the saturation level strongly depending on the 
sample size. However, also in the high-T region differ- 
ences exist: the inset to Fig.Q](b) presents L dependence 
of To for Sample I at zero magnetic field. The observed 
linearity of T) vs. hiT is consistent with the JJN model, 
and invalidates the concept of resistivity as a geometry 
independent characteristic of the TiN film. It resembles 
that of the InO K films in and strongly supports the 
idea of a collective Coulomb energy introduced in Q . Us- 
ing Eq. 12 of the data provide a single island charging 
energy E cl /k B = 0.2 K. 
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FIG. 3. (color online) Normalized I-V-curves at 28 mK and 
different sample sizes. The lines indicate power law behavior 
for a = 1 for L — 0.5, 5, and 20 /im, and a ~ 5 for 240 /im. 
Inset: R(L) at higher temperatures. 



Our central result is presented in Fig. [5] where we 
plot representative log I vs. logV traces for sample S 
(B = 0.9 T), and Sample I (B = T). In either sample 
I-Vs display switching from the highly resistive (HR)- 
to low resistive (LR) state at the threshold voltage Vhl 
(we follow convenient notations introduced in [lol]). The 
current jumps are seen over rather wide i?-range (from 
0.3 to 1.5 T for Sample S). Note that these fields are 
well below the upper critical field in TiN superconduct- 
ing films, -B C 2(0) = 2.8 T [5lj]. At higher temperatures 
T > 70 mK in both samples I-Vs are linear (ohmic) at 
low voltages but develop a characteristic steep, but con- 
tinuous increase with increasing voltage. At high volt- 
ages all I-Vs for different temperatures merge. Below 
T k 50 mK the switching voltages increase on average, 
but begin to scatter randomly over a finite temperature 
dependent V- interval (see below). 

Most striking is that both samples exhibit a power- 
law I oc V a behavior with a dramatically growing expo- 
nent a(T) in the pre-switching regime (see the inset in 
Fig. Ek). E.g. for sample S a = 2 at T = 42 mK, reach- 
ing the value a(T BKT ) = 3 at T = 37 mK and shooting 
up beyond seven at T < 32 mK. At T < 32 mK the HR 
parts of I-Vs are concealed under the noise floor and a 
finite voltage springs up at Vhl hi the LR part of the I-V 
curves. Such an evolution of the I-V curves with decreas- 
ing temperature and the corresponding T-dependence of 
a is characteristic of a charge BKT-transition which is 
expected to occur for Josephson networks in the insulat- 
ing regime, as described in the introduction 0, @ , since 
the number of free charges depends on voltage, leading 
to an intrinsically non-linear I-V characteristics. Thus 
our present findings lend strong support to the idea 0, HI 
that at very low T a crossover to a supcrinsulating state 
exists, which is of the charge BKT nature. 



We now address the size dependence of the I-V char- 
acteristics at low temperature shown in Fig. [3l For bet- 
ter comparison the / and V axis are normalized by the 
the sample size L. In absence of long range Coulomb 
effects the sheet resistance should be size independent 
and all curves in Fig. [3] should collapse. This is clearly 
not the case: the I-V characteristics for the three the 
shortest samples is linear at low V , but with a very dif- 
ferent resistivity. The size dependence of the I-V curves 
and the saturation of the sheet resistance (cf. Fig. [2Jd) 
at the lowest temperature depends is much more dra- 
matic than expected from the weak logarithmic size de- 
pendence of the activation energy. At higher V the a 
step-like feature occurs, which is most probably heating 
induced. The corresponding electric field increases with 
decreasing sample size. Only the largest sample displays 
a strongly non-linear I-V characteristic at low V with 
a very high power-law exponent a ~ 5. The inset dis- 
plays R{L) at higher temperatures T > 4 K, where the 
size dependence is reverted: R decreases with L. While 
the precise reason for this decrease is not yet clear, it 
is plausible that the patterning process, which has al- 
ready driven the sample towards the insulating side of 
the SIT, affects in particular the sample edges. A locally 
suppressed conductivity at the sample edges affects more 
strongly the smaller samples. An overall gradient of the 
local conductivity over the whole chip can be excluded, 
since the smallest and the largest sample are placed next 
to each other. 

We now discuss the properties of the sharp current 
jumps between the HR and the LR state. The latter 
had been attributed to a dielectric breakdown along a 
quasi-lD weakest dielectric strength channel 0. While a 
random distribution of Vhl is typical for dielectric break- 
down phenomena [22| . a complete description of thresh- 
old behavior requires in addition the consideration of 
heating instabilities [ill, mentioned in the introduc- 
tion, which usually accompany a dielectric breakdown. 
Such an electron overheating governs in particular the re- 
turn from the LR to the HR state upon reducing the bias 
voltage. Adopting the model considerations in Ref. 
we find that we can fit the I-V characteristics of Fig. [2] 
for T > 60 mK [24]. In contrast, at T < 60 mK the heat- 
ing model significant overestimates the range of measured 
Vhl, indicating that the heating instability is preempted 
by a dielectric breakdown at the lowest T. 

To investigate the statistical distribution of the switch- 
ing voltages Vhl at T < 60 mK, we have measured his- 
tograms of the I-V characteristics for each representa- 
tive temperature and external magnetic field. Figure 0] 
displays the corresponding histograms of Vhl for Sam- 
ple S [panels (a) and (b)]; panels (c) and (d) show the 
temperature and the magnetic field dependence of the 
maximal and minimal values of Vhl- We observe that 
i) the histograms display a sharp cut off at the maxi- 
mal observed Vhl and is asymmetric with respect to the 



4 





" maX f K HL' 

o mm{Vi 



30 40 
T [mK] 



max{F HL } JL 
g 6|-d min{C HL } 
4 



(J) 



S[T] 



FIG. 4. (color online) Histograms for the Vhl distribution: 
(a) = 0.9 T and T — 28mK (the total sampling number is 
930) and T = 55 mK (296 samplings); (b) T = 28 mK and 
B = 0.3T (776 samplings), 0.5 T (299 samplings), 1.1 T (708 
samplings), and 1.3 T (300 samplings), (c)-(d): Temperature 
(at = 0.9 T) and field (at T = 28 mK) dependences of the 
variance range of Vhl, maximal and minimal values shown by 
solid and open symbols, respectively. Horizontal lines indicate 
intervals containing 75% of the data. The solid line marked 
by Vhl corresponds to the conventional overheating model. 



most probable Vhl (B) (see the representative histogram 
for B = 0.9 T and T = 28 mK). ii) The field-dependence 
of the most probable value of Vhl(-B) is non-monotonic 
and has the values close to those observed on samples 
which are insulating at B = @. hi) The ratio of 
max{VHL(-B)}/min{VHL(S)} « 2 in magnetic fields be- 
low the field where the maximum of Uhl(-B) is achieved, 
and drops rapidly upon further field increase, iv) As tem- 
perature grows, max{VHL}/ m hi{^HL} shrinks to unity 
and the random scattering of Vhl vanishes. Note that 
the variance of Vhl disappear much faster than Vhl it- 
self. The observed jitter in the threshold voltage closely 
resembles the dielectric breakdown behavior in conven- 
tional semiconductors. A similar random switching I-V 
behavior with the characteristic current jumps and hys- 
teresis has been also observed in 2D Josephson junction 
networks in the insulating state (25- 27|. 

In conclusion the above observations reveal two char- 
acteristically different regimes within the insulating side 
of the SIT in our TiN films: i) at comparatively high tem- 
peratures 60 mK < T < 500 mK transport is linear with 
a thermally activated resistance. The activation energy 
increases with the sample size. The evolution with mag- 
netic field and size can be consistently explained within 
a Josephson network model in the insulating regime. As 
local Cooper pairing is essential for the insulating be- 
havior, we call this regime a Cooper pair insulator, ii) at 



very low temperature T < 50 mK the conduction process 
becomes intrinsically non-linear at low bias voltage. 

Hence, our low-temperature data are in accord with 
the existence of the low temperature charge Bcrezhinskii- 
Kosterlitz-Thouless phase, supporting the notion of su- 
perinsulation as dual to superconductivity. Measure- 
ments of the dielectric properties could provide indepen- 
dent and complementary evidence for this new state of 
matter. 
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